Reachability and Termination Analysis of 
Concurrent Quantum Programs 



Nengkun Yu and Mingsheng Ying 

Tsinghua University, China 
University of Technology, Sydney, Australia 
{nengkunyuOgmail . com , Mingsheng . YingOuts . edu . au} 



Abstract. We introduce a Markov chain model of concurrent quantum 
programs. This model is a quantum generalization of Hart, Sharir and 
Pnueli's probabilistic concurrent programs. Some characterizations of the 
reachable space, uniformly repeatedly reachable space and termination 
of a concurrent quantum program are derived by the analysis of their 
mathematical structures. Based on these characterizations, algorithms 
for computing the reachable space and uniformly repeatedly reachable 
space and for deciding the termination are given. 
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1 Introduction 

Research on concurrency in quantum computing started about 10 years ago, and 
it was motivated by two different requirements: 

— Verification of quantum communication protocols: Quantum communication 
systems are already commercially available from Id Quantique, MagiQ Tech- 
nologies, SmartQuantum and NEC. Their advantage over classical commu- 
nication is that security is provable based on the principles of quantum me- 
chanics. As is well known, it is very difficult to guarantee correctness of even 
classical communication protocols in the stage of design. Thus, numerous 
techniques for verifying classical communication protocols have been devel- 
oped. Human intuition is much better adapted to the classical world than 
the quantum world. This will make quantum protocol designers to commit 
many more faults than classical protocol designers. So, it is even more criti- 
cal to develop formal methods for verification of quantum protocols (see for 
example [TU], [H], [1]). Concurrency is a feature that must be encompassed 
into the formal models of quantum communication systems. 

— Programming for distributed quantum computing: A major reason for dis- 
tributed quantum computing, different from the classical case, comes from 
the extreme difficulty of the physical implementation of functional quan- 
tum computers (see for example [I], |21|). Despite convincing laboratory 
demonstrations of quantum computing devices, it is beyond the ability of 



the current physical technology to scale them. Thus, a natural idea is to use 
the physical resources of two or more small capacity quantum computers to 
simulate a large capacity quantum computer. In fact, various experiments in 
the physical implementation of distributed quantum computing have been 
frequently reported in recent years. Concurrency naturally arises in the stud- 
ies of programming for distributed quantum computing. 

The majority of work on concurrency in quantum computing is based on 
process algebras E], E], [E], 0, E], 0, E2, 0, [3]. This paper introduces a 
new model of concurrent quantum programs in terms of quantum Markov chains. 
This model is indeed a quantum extension of Hart, Sharir and Pnueli's model of 
probabilistic concurrent programs |12| . |19j . Specifically, a concurrent quantum 
program consists of a finite set of processes. These processes share a state Hilbert 
space, and each of them is seen as a quantum Markov chain on the state space. 
The behaviour of each processes is described by a super-operator. This descrip- 
tion of a single process follows Selinger, D'Hont and Panangaden's pioneering 
works |18j . [S] on sequential quantum programs where the denotational semantics 
of a quantum program is given as a super-operator. The super-operator descrip- 
tion of sequential quantum programs was also adopted in one of the authors' 
work on quantum Floyd-Hoare logic [30] • Similar to the classical and probabilis- 
tic cases |12| . an execution path of a concurrent quantum program is defined 
to be an infinite sequence of the labels of their processes, and a certain fairness 
condition is imposed on an execution path to guarantee that all the processes 
fairly participate in a computation. 

Reachability and termination are two of the central problems in program 
analysis and verification. The aim of this paper is to develop algorithms that 
compute the reachable states and decide the termination, respectively, of a con- 
current quantum program. To this end, we need to overcome two major diffi- 
culties, which are peculiar to the quantum setting and would not arise in the 
classical case: 

— The state Hilbert space of a quantum program is a continuum and thus 
doomed-to-be infinite even when its dimension is finite. So, a brute-force 
search is totally ineffective although it may works well to solve a correspond- 
ing problem for a classical program. We circumvent the infinity problem of 
the state space by finding a finite characterization for reachability and termi- 
nation of a quantum program through a careful analysis of the mathematical 
structure underlying them. 

— The super-operators used to describe the behaviour of the processes are op- 
erators on the space of linear operators on the state space, and they are 
very hard to directly manipulate. In particular, algorithms for computing 
super-operators are lacking. We adopt a kind of matrix representation for 
super-operators that allows us to conduct reachability and termination anal- 
ysis of quantum programs by efficient matrix algorithms. 

The paper is organized as follows. For convenience of the reader we briefly 
recall some basic notions from quantum theory and fix the notations in Sec. [5] but 



we refer to |17) for more details. A Markov chain model of concurrent quantum 
programs is defined in Sec. [3J where we also give a running example of quantum 
walks. In Sec. SI we present a characterization for reachable space and one for 
uniformly repeatedly reachable space of a quantum program, and develop two 
algorithms to compute them. A characterization of termination of a quantum 
program with fair execution paths and an algorithm for deciding it are given in 
Sec. [5] It should be pointed out that termination decision in Sec. [5] is based on 
reachability analysis in Sec. |U A brief conclusion is drawn in Sec. |51 

2 Preliminaries and Notations 

2.1 Hilbert Spaces 

The state space of a quantum system is a Hilbert space. In this paper, we only 
consider a finite dimensional Hilbert space "H, which is a complex vector space 
equipped with an inner product (-|-). A pure state of a quantum system is rep- 
resented by a unit vector, i.e., a vector \ip) with (ip\ip) = 1. Two vectors \<p), |i/>) 
in H are orthogonal, written \<p) _L if their inner product is 0. A basis of % 
is orthonormal if its elements are mutually orthogonal, unit vectors. The trace 
of a linear operator A on T~L is defined to be tr(A) — '^2, i (i\A\i) , where {\i}} is an 
orthonormal basis of T-L. For a subset V of H, the subspace spanF spanned by V 
consists of all linear combinations of vectors in V. For any subspace X of T-L, its 
orthocomplement is the subspace X 1 - = {\<p) 6 H : \<p) _L \ip) for all \ip) € X}. 
The join of a family {Xi} of subspaces is \/ i Xi = span(ljj Xi). In particular, we 
write X V Y for the join of two subspaces X and Y. A linear operator P is called 
the projection onto a subspace X if P\ip) = for all \ip) G X and P\ip) = 
for all \ip) £ .X" . We write Px for the projection onto X . 

A mixed state of a quantum system is represented by a density operator. 
A linear operator p on % is called a density operator (resp. partial density 
operator) if p is positive-semidefinite in the sense that {(j)\p\<f>) > for all \<f>), 
and tr(p) = 1 (resp. tr(p) < 1). For any statistical ensemble {(pi, IV'i))} °f 
pure quantum states with pi > for all i and ^2tPi = 1, p = ^iPilip^i^il i s 
a density operator. Conversely, each density operator can be generated by an 
ensemble of pure states in this way. In particular, we write ij) for the density 
operator |V>) generated by a single pure states The support of a partial 
density operator p. written supp(p), is the space spanned by its eigenvectors 
with nonzero eigenvalues. 

Lemma 1 For any p > and partial density operators p,cr, we have: (1) 
supp(pp) = supp(p); (2) supp(p) C supp(p + a); (3) supp(p + a) = supp(p) V 
supp(a). 

2.2 Super-Operators 

A super-operator is a mathematical formalism used to describe a broad class 
of transformations that a quantum system can undergo. A super-operator on 



H is a linear operator £ from the space of linear operators on H into itself, 
satisfying (1) tr[£(p)] < tr(p) for any p; (2) Complete positivity(CP): for any 
extra Hilbert space %k. (T^ ®£){A) is positive provided A is a positive operator 
on "H/c ®%, where Ik is the identity operation on H.^. Furthermore, if tr[£(p)\ = 
tr(p) for any p, then £ is said to be trace-preserving. Each super-operator £ 
enjoys the Kraus representation: there exists a set of operators {Ei} satisfying 
(1) £{p) = J^i EipEj for all density operators p; (2) £\ E\E l < I, with equality 
for trace-preserving £, where I is the identity operator. In this case, we write £ = 
J2i Ei-E\. The image of subspace X of H under £ is £ (X) = Vi-iMeX supp(£ '(ip)) , 
and the pre-image of X under £ is £ ~ 1 (X) = {\ip) E H : supp(£ (ip)) C X}. 

Lemma 2 (1) supp(p) C supp(cr) supp(£(p)) C supp(£(cr)) ; and supp(p) = 
supp(cr) => supp(£(p)) = supp(£(tx)). 

(2) su PP (£(p)) C supp((£ + r)(p)). (3) £(X) = S vl PP (£(P x )). 

(4) X <ZY ^£(X) Q£(Y). (5) £(X) C (£+F)(X). 

(6) If£ = J2i Et-Ej, then £~\X) = [supp(£ *(P x x))] x , where £* = £. E\ ■ 
Ei is the (Schrddinger-Heisenberg) dual of £ . 

2.3 Matrix Representation of Super-Operator 

The matrix representation of a super-operator is usually easier in [53] to manip- 
ulate than the super-operator itself. If £ — ^ Ei ■ E\ and dim'H = d, then the 
matrix representation of £ is the d 2 x d 2 matrix M = J^i^i ® E* , where A* 
stands for the conjugate of matrix A, i.e., A* = (a*j) with a*j being the conju- 
gate of complex number , whenever A — (a,ij ) . According to |24j , we have the 
following 

Lemma 3 (1) The modulus of any eigenvalue of M is less or equal to 1. 

(2) We write \<&) = . \ for the (unnormalized) maximally entangled state 
inT-L^T-L, where is an orthonormal basis ofH. Then for any d x d matrix 

A, we have {£ (A) <g> 7)|#) = M(A®I)\$). 

2.4 Quantum Measurements 

A quantum measurement is described by a collection {M m } of operators, where 
the indexes m refer to the measurement outcomes. It is required that the mea- 
surement operators satisfy the completeness equation J2m M^Mm = In- If the 
system is in state p, then the probability that measurement result m occurs is 
given by p(m) = tr(M} n M m p), and the state of the system after the measurement 

p(m) 

3 A Model of Concurrent Quantum Programs 

Our model is a quantum extension of Hart, Sharir and Pnueli's probabilistic 
concurrent programs |12) . A concurrent quantum program consists of a finite set 



K = {1,2, ■ • • , to} of quantum processes, and these processes have a common 
state space, which is assumed to be a <i-dimensional Hilbert space T-L. With each 
k € K we associate a trace-preserving super-operator £)., describing a single 
atomic action or evolution of process k. Also, we assume a termination condi- 
tion for the program. At the end of each execution step, we check whether this 
condition is satisfied or not. The termination condition is modeled by a yes-no 
measurement {Mo, Mi}: if the measurement outcome is 0, then the program ter- 
minates, and we can imagine the program state falls into a terminal (absorbing) 
space and it remains there forever; otherwise, the program will enter the next 
step and continues to perform a quantum operation chosen from K. 

Definition 1 A concurrent quantum program defined on a d- dimensional Hilbert 
space % is a pair V — {{£k ■ k 6 K}, {Mo, Mi}), where: 

1. £k is a super- operator on H for each k g K ; 

2. {Mo, Mi} is a measurement on T-L as the termination test. 

Any finite string s\S2 ■ • • s m or infinite string s\S2 ■ • • Sj • • • of elements of K 
is called a execution path of the program. Thus, the sets of finite and infinite 
execution paths of program V are 

S = K u — {siS2 ■ ■ ■ Si ■ ■ ■ : Si e K for every i > 1}, 
Sfin = K* = {siS2 ■ ■ ■ s m : m > and s, € K for all 1 < i < to}, 

respectively. A subset of S is usually called a schedule. 

For simplicity of presentation, we introduce the notation J-^ for any k G K 
which stands for the super-operator defined by Tk(p) = EkiM\pM\) for all den- 
sity operators p. Assume the initial state is pq. The execution of the program 
under path s = S1S2 ■ ■ ■ Sfc • • • € S can be described as follows. At the first step, 
we perform the termination measurement {Mo, Mi} on the initial state po- The 
probability that the program terminates; that is, the measurement outcome is 
0, is It[MqPqMq]. On the other hand, the probability that the program does 
not terminate; that is, the measurement outcome is 1, is p\ = ti[AlipoM\], and 
the program state after the outcome 1 is obtained is p\ — M\poM\/p\. We 
adopt Selinger's normalization convention |18| to encode probability and density 
operator into a partial density operator p\p\ — Mip M\. Then this (partial) 
state is transformed by the quantum operation £ Sl to £ Sl (MipoMi) = J 7 Sl (po). 
The program continues its computation step by step according to the path s. 
In general, the (n + l)th step is executed upon the partial density operator 
PnPn = ^ o • • • o J r S2 o J- Sl (po), where p s n is the probability that the pro- 
gram does not terminate at the nth step, and p s n is the program state after 
the termination measurement is performed and outcome 1 is reported at the 
nth step. For simplicity, let Tt denote the super-operator J- Sn o • • • o T s , 2 o T S1 
for string / = s\Sf ■ ■ s n . Thus, p s n p s n = J~ s [n](Po), where s[n] is used to de- 
note the head s\S2 ■ • • s n for any s = s\S2 ■ • • s n • • • € S. The probability that 
the program terminates in the (n + l)th step is then tr(M (.F s [ n ] (pq))Mq), and 
the probability that the program does not terminate in the (n + l)th step is 
p« + i = tr(M 1 (J- sW (p ))M 1 t ). 



3.1 Fairness 



To guarantee that all the processes in a concurrent program can fairly participate 
in a computation, a certain fairness condition on its execution paths is needed. 

Definition 2 An infinite execution path s — s\S2---Si... G S is fair if each pro- 
cess appears infinitely often in s; that is, for each k G K, there are infinitely 
many i>l such that s, — k. 

We write F — {s : s € S is fair} for the schedule of all fair execution paths. 

Definition 3 A finite execution path a = s\Si • • • s n € S/,„ is called a fair piece 
if each process appears during a; that is, for each k e K, there exists i < n such 
that Si — k. 

Ffi n is used to denote the set of all fair pieces: Ff in = {a : a e Sji n is a fair piece} 
It is obvious that F = Fj in ; in other words, every fair infinite execution path 
seF can be divided into an infinite sequence of fair pieces: s = /1/2 ■ ■ ■ fk ■ ■ ■ j 
where fi e Ff in for each i > 0. The fairness defined above can be general- 
ized by introducing the notion of fairness index, which measures the occurrence 
frequency of every process in an infinite execution path. 

Definition 4 For any infinite execution path s G F, its fairness index f(s) is 
the minimum, over all processes, of the lower limit of the occurrence frequency 
of the processes in s; that is, 



where s(n, fc) is the number of occurrences of k in s[n\. 

For any S > 0, we write Fg for the set of infinite execution paths whose fairness 
index is greater than 5: F$ — {s : s e S and f(s) > 5}. Intuitively, within 
an infinite execution path in Fg, each process will be woken up with frequency 
greater than S. It is clear that F C F. 

3.2 Running Example 

We consider two quantum walks on a circle C3 = (V, E) with vertices V = 
{0,1,2} and edges E = {(0, 1), (1, 2), (2, 0)}. The first quantum walk Wi = 
({Wi}, {M , Mi}) is given as follows: 

— The state space is the 3— dimensional Hilbert space with computational basis 



— The initial state is |0); this means that the walk starts at the vertex 0; 

— A single step of the walk is defined by the unitary operator: 



where w = e 2 W3_ Intuitively, the probabilities of walking to the left and to 
the right are both 1/3, and there is also a probability 1/3 of not walking. 



f(s) = min lim inf 

fee K t^oo n>t 



s(n, k) 



n 



{\i)\i <E V}; 



Wi = — ■= 1 w w 2 




— The termination measurement {Mo, Mi} is defined by 



M = |2)(2|, M 1 =J 3 -|2)<2|, 
where ^3 is the 3x3 unit matrix. 

The second walk W2 = {{W2}, {Mo, Mi}) is similar to the first one, but its single 
step is described by unitary operator 

1 I 1 1 M 

W 2 = -p 1 w 2 w . 
V 3 \lww 2 J 

Then we can put these two quantum walks together to form a concurrent pro- 
gram V = ({Wi, W2}, {Mo, Mi}). For example, the execution of this concurrent 
program according to unfair path 1" ^ F is equivalent to a sequential program 
({Wi}, {Pq,Pi}); and the execution of V according to fair path (12)" G F is as 
follows: we perform the termination measurement {Mo, Mi} on the initial state 
po, then the nonterminating part of the program state is transformed by the 
super-operator Ui = W\ ■ W\, followed by the termination measurement, and 
then the application of the super-operator U2 — W2 ■ W\ , and this procedure is 
repeated infinitely many times. 



4 Reachability 

Reachability is at the centre of program analysis. A state is reachable if some 
finite execution starting in the initial state ends in it. What concerns us in the 
quantum case is the subspace of T-L spanned by reachable states. 

Definition 5 The reachable space of program V = {{£k '■ k £ K}, {Mo, Mi}) 
starting in the initial state po is 

H R = \f supp J" sb ](po) = \/ supp Ff{p ). 

seS,j>o fes fin 

We have the following closed form characterization of the reachable space. 

Theorem 1 Hjj = supp(^^ =Q 1 J-^po)), where d = dim"H is the dimension of 
U, and T = J^keK^k- 

Proof: We write X for the right-hand side. From Lemma [TJ we see that X = 
V{suppJ r /( / 9 ) : / G Sfi n , l/l < d}, where |/| denotes the length of string 
/. According to the definition of reachable space, we know that X C % R . To 
prove the inverse part X D Hr, for each n > 0, we define subspace Y n as 
follows: Y n := supp(^" =0 F l (po))- Due to Lemma [U we know that Yo C Yi C 
■ ■ • C Y n C • ■ ■ . Suppose r is the smallest integer satisfying Y r — Y r+ i. We 
observe that Y n+ i = supp(p + ^(Py^)) for all n > 0. Then it follows that 
Yn = Yr for all n > r. On the other hand, we have Yq £ Y\ C • • ■ C Y r . 



So, < do < d\ < ■ ■ ■ < d r < d, where do is the rank of po, and di is the 
dimension of subspace Yj, for 1 < i < r. Therefore, we have r < d — 1 and 
Yd-\ = Y r D Y n for all n. Finally, for any / £ Sfi n , it follows from Lemma [2] 
that supp(J/(p )) C supp^^Oo)) C Y\ f \ C Y d ^ = X. Thus, H R CI ■ 
Now we are able to present an algorithm computing reachable subspace using 
matrix representation of super-operators. We define Q — ^2 kGK !Fk/\K\. 



Algorithm 1: Computing reachable space 
input : An input state po, and the matrix representation G of Q 
output: An orthonormal basis B of T-Lr. 

(* \$) — ^2 . \ja3b) is the unnormalized maximally entangled state in H ® H *) 
for j = 1 : d do 

I \Vi) <~ 0'b|x); 
end 

set of states B 0; 
integer / <s— 0; 
for j = 1 : d do 

k> «- l%> - Ei=i( & fc|yj>l & fc>; 

if |z} / then 

I -e 2 + 1; 

|6j) <- ^/^/(ili); 
B^BU {»}; 
end 
end 
return 



Theorem 2 Algorithm 1 computes the reachable space in time 0(<i 4 - 7454 ) , 
where d — dimH. 

Proof: It follows from Lemma[3Jl) that I—G/2 is invertible, and X^o^/^)* = 
(/- G/2)-\ We write p = J2Zo SHpa)/2 l , and have 

(p®/)|#> = ^(G/2)>o® /)!<?>), 

i=0 

and the existence of p immediately follows from Lemma [3] We further see that 
|x) = (p® I)\$) = ~Y^j p\3a)\jb) and = p\jA)- Note that _B is obtained from 
{|%)} by the Gram-Schmidt procedure. So, supp(p) = span{p|j)} = spamB. It 
is clear that Hr = supp(^^T 1 J"*(po)) Q su PP(p)- Therefore, Hr = supp(p) = 
spanB, and the algorithm is correct. 

The complexity comes from three the following parts: (1) it costs 0(d 2 * 2 ' 3727 ) 
to compute (I — G/2) -1 by using Coppersmith- Winograd algorithm [2J; (2) it 



requires 0(d A ) to obtain \x) from (I — G/2) -1 ; (3) the Gram-Schmidt orthonor- 
malization is in time 0(d 3 ). So, the time complexity is C(d 4 7454 ) in total. ■ 
An advantage of Algorithm 1 is that we can store (I — G/2) -1 . Then for any 
input state po, we only need 0(d ) to compute the space reachable from po- 

Definition 6 The uniformly repeatedly reachable space of program V = ({£ k ■ 
k G K}, {Mo, -Mi}) starting in the initial state po is 

Hurr = P| \J supp F s \j](po) = P| V"t su P p Ffipo) ■ f G 5 /m , \f\ > n}. 

n>OseSj>n n>0 

The uniformly repeatedly reachable space enjoys the following closed form, 
Theorem 3 Hurr — supp(J^ = ~ d 1 ^{po)), where d = dimH, and J 7 = J2kei< -^k 

Proof: For each n > 0, we define subspace Z n as follows: Z n := V ' - >n supp ^(po). 
It is obvious that 3 2i 3 ■ ■ • 3 Z„ D • • • . Suppose r is the smallest integer 
satisfying Z r = Z r+ i. By noting that Z n+ i = supp(J r (P2 n )), we can show that 
Z n = Z r for all n > r. On the other hand, we have Z$ ^> Zi D ■ ■ ■ D Z r . So, 
do > d\ > ■ ■ ■ > d r > 0, and di is the dimension of subspace Zi for < i < r. 
Therefore, we have r < do < d and Zd = Z r . Therefore, Hurr = HnX) ^ n = 
It is obvious that Zd is the reachable space starting in state T {po)- Using The- 
oremQ] we obtain Z d = supp(X^=o F l (F d (po))) = suppQ^Lj 1 F(po))- ■ 
We can give an algorithm computing the uniformly repeatedly reachable space 
by combining the above theorem and matrix representation of super-operators. 



Algorithm 2: Compute uniformly repeatedly reachable space 
input : An input state po, and the matrix representation G of Q 
output: An orthonormal basis Burr of Hurr- 
\x)^G d {I-G/2)- x {po®I)\$); 

(* \4>) = ^2 j \ja3b) is the unnormalized maximally entangled state in H <£> Tt *) 
for j = 1 : d do 

I \Vi) <~ 0'b|x); 
end 

set of states Burr 0; 
integer / <s— 0; 
for j = 1 : d do 

l*)Hvi>-£U< 6 *lw>lfc>; 

if \z) ± then 

J <- 2 + 1; 

Burr <- Bras U{|(n)}; 
end 
end 
return 



Theorem 4 Algorithm 2 computes the uniformly repeatedly reachable space in 
time 0(d 4 7454 log d), where d — dim'H. 

Proof: This theorem is a corollary of Theorem [21 Here, logd in the complexity 
comes from computing M d using the method of exponentiation by squaring. 

5 Termination 

Another important problem concerning the behaviour of a program is its termi- 
nation. 

Definition 7 Let the program V = ({£ k ■ k € K},{M ,Mx}). Then V with 
input po terminates for execution path s G S if p s n p s n = J~ s [n] = f or some 
positive integer n. 

Definition 8 1. If a program V with input po terminates for all s G A, then 

we say that it terminates in schedule A. 
2. If there is a positive integer n such that p s n p s n — for all s G A, then it is 
said that the program V with input po uniformly terminates in schedule A. 

We hrst prove the equivalence between termination and uniform termination. 
Of course, this equivalence comes from finiteness of the dimension of the state 
space. 

Theorem 5 The program V = {{£k '■ k G K }, {Mo, Mi}) with initial state po 
terminates in the biggest schedule S — K u if and only if it uniformly terminates 
in schedule S . 

Proof. The "if" part is obvious. We prove the "only if" part in two steps: 

(1) We consider the case of \K\ = 1, where : k G K} is a singleton {£}. 
Now the program is indeed a sequential program, and it is a quantum loop |23| . 
We write J-(p) = £{MipM{) for all p. What we need to prove is that if V 
terminates, i.e., T n \p§) — for some n, then it terminates within d steps, i.e., 
T d (po) — 0. If po is a pure state then we define the termination sets as 
follows: X n := {\ip) : J- n (t^) = 0} for each integer n > 0. 

(1.1) If \<p),\x) G X n , then T n (ip + X ) = 0, which leads to a\<p) + (3\x) G X n 
for any a, f3 G C. Thus X n is a subspace of H. 

(1.2) Since F n {^) = T n+1 (ijj) = 0, it holds that that X n C X n+1 for 
any n > 0. So, we have the inclusion relation X\ C X2 C • • • C X n C • • • . 

Now suppose t is the smallest integer satisfying X t — X t+ \. Invoking Lemma 
H we obtain that supp^*^/)) = = X^ +l = supp^'+^I)), where 
denotes the (Schrodinger-Heisenberg) dual of We have supp(J r *"(/)) = 

supp (J 7 * (7)), which leads to X n = X t for all n > t. Now, it holds that X\ C 
X2 C • • • C Xt = X t +\ — X t+ 2 — • • ■ • This implies dx < d% ■■ ■ < dt, where di is 
the dimension of subspace Xi. Thus, t < d. If J 711 ^) = 0, then \ip) G X n C Xj, 
and F d {iJ}) = 0. 



In general, if po is a mixed input state po = '%2Pi\' t l ) i){' , Pi\ with all pi > 0, and 
F n (p ) = 0, then F n (ipi) = for all i. Therefore, = for all i, and it 

follows immediately that T (po) = 0. 

(2) For the general case of \K\ > 2, we assume that V starting in po terminates 
in S, i.e., for any s £ S, there exists an integer n s such that J- s [ ns ](po) = 0. Our 
purpose is to show that there exists an integer n such that J- s [ n ](po) = f° r 
all s £ S. Indeed, we can choose n — d. We do this by refutation. Assume that 
■F s [d]{po) 7^ for some s £ S. We are going to construct an execution path 
s £ S such that ^F s [n](po) ^ for any n > 0. Let T = J2k£K J~k- Then the 
assumption means that there exist / £ K d such that Ff(po) ^ 0, and it follows 
that F d (po) 7^ 0. Now we consider the loop program {{J 7 }, {Mo, Mi}) with initial 
state po. Applying (1) to it, we obtain F 2d (po) ^ 0. Then there exist g\, hi £ K d 
such that FhxiFgAPo)) = Fgihx(po) ^ 0, and F d {F gi (po)) 7^ 0. Applying (1) 
again leads to F 2d (F gi (po)) 7^ 0, which means that there exist hi, 52 £ K 2d such 

that J r h 2 (J r g 1 g 2 (Po)) = FgihaiFgiiPo)) ^ °- Tllus : we nave ^(-^flafli (po)) ^ 0. 

Repeating this procedure, we can find an infinite sequence gx, 52, •■• £ K d . Put 
s = 9l92— £ S. Then it holds that T s [kd](Po) 7^ for any integer k. Thus, we 
have T s [n] (p) 7^ for all n. ■ 

Now we are ready to consider termination under fairness. Of course, any per- 
mutation of K is a fair piece. We write Pk for the set of permutations of K. For 
a = s\Si ■ • • s m £ Pk, a finite execution path of the form siCi^o^ ■ ■ ■ 0Vn-l s m 
is called an expansion of a. Obviously, for any a £ Pk, all of its expansions are 
in Ffi n . We will use a special class of fair pieces generated by permutations: 

n = {siU\S2Cf2 ■ ■ ■ Grn-\Srn. '■ S1S2 ■ ■ ■ s m £ Pk and \<Ji\ < d for every 1 < i < m}, 

where d is the dimension of the Hilbert space H of program states. It is easy to 
see that II C Ffi n . 

Theorem 6 A program V = ({St '■ k £ i^},{Mo,Mi}) with initial state po 
terminates in the fair schedule F if and only if it terminates in the schedule 77". 

Proof. The "only if" part is clear because 77" C F. To prove the "if" part, assume 
V terminates in the schedule 77". We proceed in four steps: 

(1) Since 77 is a finite set, we can construct a new program V' = ({J~f ■ f £ 
77}, {0,7}). (We should point out that Fj is usually not trace-preserving, and 
thus T" is indeed not a program in the sense of Definition \T\ However, this does 
not matter for the following arguments.) It is easy to see that the termination 
of V with po in schedule 77" implies the termination of V with po in 77". Note 
that 77" is the biggest schedule in V' , although it is not the biggest schedule 
in V . So, we can apply Theorem [S] to V' and assert that (J2fen-^f) d (Po) = 
That is equivalent to 




(1) 



fen 



(2) For each a G Pk, we set A„ = {a 1 6 77 : a' is an expansion of <r}. Then 
(J A a = 77. Moreover, we write Q a = J^feA J~f f° r ever y a £ Pk- It is worth 

noting that ~^2 a( zp K Q a = S/gtt ^/ ^ s n0 ^ true in general because it is possible 
that A ai n A CT2 ^ for different o\ and 02 ■ But by Lemma [T]l and 3 we have 

SUpp[( ^ &<r)(f>0)] = SUpp[(J^ Ff)(po)], 

aeP K fen 
and furthermore, it follows from Eq. (TT]) that 

supp[( Sa) d ( P o)} = supp[(£ J» d (po)] = {0}. (2) 
aeP K fen 

(3) For each fair piece a' G Ffi n , and for any p, we can write a' = s±fiS2 ■ ■ ■ 
^m-l/m-l^m for some cr = sis 2 ---s m e Pk, and /i,...,/ m _i € S/ m . Fur- 
thermore, we write £/ = 5Zi=(jG^fcLi -^"fc) 1 - First, a routine calculation leads to 
So-o = ?s m ° ° Js,,-! • • • -7 r 52 ° 5 ° ^"si- Second, it follows from Theorem [T] that 
for each 1 < i < m — 1, and for any p, supp(J r / i (/>)) C supp(^( j o)). Repeatedly 
applying this inclusion together with Lemma [2] 1 we obtain 



SUpp(Jv/(p)) = SUpp[(J" Sm o J 7 /^ o J r Sm _ 1 o • • • o J" S2 o T h o J" Sl )(p)] 
C supp[(J" Sm o g o J" Sm _j o.-oj 8j ogo T Sl )(p)] 

= supp(£ CT0 (p)) C supp(J^ -7 r <r)(p)- 

cr£77 



(3) 



(4) Now we are able to complete the proof by showing that for any fair 
execution path s G F, s has an initial segment t such that Tt(po) = 0. In fact, s 
can be written as an infinite sequence of fair piece, i.e., s = o-[a' 2 • • • , where each 
a'i is a fair piece. We take t to be the initial segment of s containing the first d 
fair pieces, i.e., t = o\o<2, ■ ■ ■ ad- Repeatedly applying Eq. ((3]) and Lemma [2] 1 we 
obtain 

suppJ" t (p ) = supp[(J^ o-oj.,0 -TvjXpo)] 

Csupp[(£ F p ) d (p)} = {0}. 
pen 

Thus, T t {p) = 0. ■ 

The above theorem can be slightly strengthened by employing the notion of 
fairness index in Definition [U First, we have: 

Lemma 4 77" C Fx . 

Proof. For any s = <j\Oi • • • G 77" with er,; G 77, we know that <7j(|<7i|, k) > 1 for 
any k G K and < md, where cr^ ( | cr^ | , k) is the number of occurrences of k 
in <7j. Then the occurrence frequency f(s)> which means that 77" C F i . 
On the other hand, we choose an arbitrary s G F i . Then l md s G F i but 
l md s <^ 77". ■ 



Actually, what we proved in Theorem [5] is that for any two schedules A, B 
between 27" and F, i.e., 77" C A, B C F, a program terminates in schedule A 
if and only if it terminates in schedule B . Combining Theorem [5] and Lemma 0] 
yields: 

Corollary 1 For any < S, e < -^g, a program terminates in schedule F$ if and 
only if it terminates in schedule F e . M 

Now an algorithm checking termination of a concurrent quantum program 
can be developed based on Theorem [SJ 



Algorithm 3: Decide termination of a concurrent quantum program 
input : An input state po, and the matrix representation of each Ti i.e, Ni 
output: b.(If the program terminates under F, b = 0; Otherwise, b — 1.) 
N <- 0; 

for k = 1 : m do 

| N<-Ni+N; 
end 
G^I; 

for fe = 1 : d — 1 do 
| G 4— I + NG; 
end 

(*Compute the matrix representation of Q*) 
M <~ 0; 
Generate Pk\ 

for p = pip 2 ■ ■ -p m G Pk do 
L<-N P1 ; 
for I = 2 : m do 
| L 4- N Pl GL; 
end 

(^Compute the matrix representation of J- p *) 
M<-M + L; 
end 

(*Compute the matrix representation of ^2 pe p K F p *) 

\x) ^Af d (po® /)|<£); 
if \x) / then 
I b<-l; 
end 

if [a;) = then 

I fe^O; 
end 

return b 



Theorem 7 Algorithm 3 decides termination of a concurrent quantum program 
in time 0(m m d 4 ' 7454 ), where m is the number of the processes, and d = dim%. 



Proof: In the algorithm, we use the for loop to compute the matrix rep- 
resentation G of Q = Y^i=o(Y^k=i -^k) 1 - Then the matrix representation of 
T a = J- Sl o Q o ■ ■ • Q o T Sm (•) is obtained for any a — s\S2 ■ ■ ■ s m G Pk . 
All T a s are added up to M. Then M becomes the matrix representation of 
Sa-ePjf Consequently, we can apply TheoremlHlto assert that this algorithm 
outputs if the program terminates in the fair schedule F; otherwise, 1. 

To analyse its complexity, the algorithm can be divided into three steps: (1) 
Computing G costs 0(m + d d 2 * 2 - 3727 ) = 0(m + d 5 - 7454 ); (2) Computing M costs 
m!*2m*0(d 2 * 2 - 3727 ) = 0(m m d 4 - 7454 ); (3) Computing \x) costs 0{d A - 7454 log d). 
So the total cost is 0((m m + d)d 4 - 7454 ). U 

6 Conclusion 

In this paper, we studied two of the central problems, namely, reachability and 
termination for concurrent quantum programs. A concurrent quantum program 
is modeled by a family of quantum Markov chains sharing a state Hilbert space 
and a termination measurement, with each chain standing for a participating 
process. This model extends Hart, Sharir and Pnueli's model of probabilistic 
concurrent programs [T2] to the quantum setting. We show that the reachable 
space and the uniformly repeatedly reachable space of a concurrent quantum 
program can be computed and its termination can be decided in time 0(d 4 7454 ), 
0(<i 4 ' 7454 log d), 0((m m + d)cZ 4 ' 7454 ), respectively, where m is the number of 
participating processes, and d is the dimension of state space. 

For further studies, an obvious problem is: how to improve the above al- 
gorithm complexities? In this paper, reachability and termination of quantum 
programs were defined in a way where probabilities are abstracted out; that is, 
only reachability and termination with certainty are considered. A more delicate, 
probability analysis of the reachability and termination is also an interesting 
open problem. The algorithms for computing the reachable space and checking 
termination of a quantum program presented in this paper are all algorithms for 
classical computers. So, another interesting problem is to find efficient quantum 
algorithms for reachability and termination analysis of a quantum program. 
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